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Abstract: We present first neutrino induced events observed with a deep underwater neutrino 
telescope. Data from 70 days effective life time of the BAIKAL prototype telescope NT-96 have 
been analyzed with two different methods. With the standard track reconstruction method, 9 
clear upward muon candidates have been identified, in good agreement with 8.7 events expected 
from Monte Carlo calculations for atmospheric neutrinos. The second analysis is tailored to 
muons coming from close to the opposite zenith. It yields 4 events, compared to 3.5 from Monte 
Carlo expectations. From this we derive a 90% upper flux limit of 1.1 • 10 -13 cm -2 sec -1 for 
muons in excess of those expected from atmospheric neutrinos with zenith angle > 150 degrees 
and energy > 10 GeV. 
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1 Introduction 



The ultimate goal of large underwater neutrino telescopes is the identification of extrater- 
restrial neutrinos of high energies. The omnipresent background with respect to these 
neutrinos are neutrinos generated by cosmic ray interactions in the atmosphere of the 
Earth. While being a background with respect to extraterrestrial neutrinos, atmospheric 
neutrinos can be used as a standard signal to test and calibrate underwater neutrino de- 
tectors. Therefore, a primary challenge for these detectors is the identification of upward 
muons generated in interactions of atmospheric neutrinos. Taking into account that the 
flux of downward muons at 1 km depth is about 6 orders of magnitude larger than the 
flux of upward muons from atmospheric neutrinos, this task is far from being trivial JTJ. 

Apart from being a test tool for underwater telescopes, atmospheric neutrinos may indi- 
cate physics beyond the standard model of particle physics. Actually, from recent analy- 
ses of data from underground experiments, evidence is found for oscillations of neutrinos 
[0, ||, ||, £§]. The search for such feeble effects, however, requests a detailed understanding 
of the detectors, high statistics and a relatively low energy threshold. 

The first part of this paper describes the method for full spatial reconstruction of muon 
tracks with the four-string array NT-96 in Lake Baikal, and the criteria to reject fake 
events from misreconstructed downward muons. We present results from the analysis of 
70 days effective lifetime of NT-96, taken in 1996. 

In contrast to the standard reconstruction strategy, which supposes >6 hits at >3 of the 
vertical strings, the second analysis is performed also for events with hits on less than 
3 strings. This considerably increases the effective area in vertical direction. Instead of 
beginning with a reconstruction and then applying quality cuts, we start with cuts effec- 
tively rejecting all events with the exception of nearly vertically upward moving muons. 
Only for events surviving these cuts, a fit of the zenith angle is performed. An excess of 
such muons over the expectation value for muons from atmospheric neutrinos could indi- 
cate neutrino production by the annihilation of neutralinos - the favored super-symmetric 
candidate for cold dark matter - in the center of the Earth. By a method of that kind, 
two first neutrino events have been identified in NT-36, the very first small BAIKAL 
prototype array 0. In the second part of this paper, we apply this strategy to NT-96 
data and derive an upper limit on the excess muons from the center of the Earth. 
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2 The detector 



The Baikal Neutrino Telescope QTJ, H is operated in the Siberian Lake Baikal, at a depth 
of 1.1km. In April 1998, the NT-200 array consisting of 192 optical modules (OMs) 
at 8 strings was deployed (see fig.0). The OMs are grouped in pairs along the strings. 
They contain 37-cm QUASAR phototubes |J. The two tubes of a pair are switched in 
coincidence in order to suppress background from bioluminescence and dark noise. A pair 
defines a channel. 




72 m 



Figure 1: Sketch of the Baikal Neutrino Telescope NT-200. Top right, the 1996 stage NT-96 is 
shown. 

A muon trigger is formed by the requirement of > N hit channels within 500 nsec. N 
is typically set to the value of 3 or 4. For such events, amplitude and time of all fired 
channels are digitized and sent to shore. The event record includes all hits within a time 
window of -1.0 /isec to +0.8 /isec with respect to the muon trigger signal. 

In April 1993, the first small prototype array, the detector NT-36 with 36 OMs at 3 
short strings, was launched and took data up to March 1995. Results obtained with the 
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1993 version of NT-36 as well as a detailed description of detector, site and deployment 
procedure have been published in A 72-OM array, NT- 72, run in 1995-96. In 1996, it 
was replaced by the four-string array NT-96, which was replaced a year later by a 144-OM 
detector. Since April 6, 1998, NT-200 with 192 OMs is taking data. 

NT-96 (48 channels) is shown at top right of fig.[l]. The OMs are fixed at four strings each 
72 m long. The strings are arranged at the edges of a trapezoid with side lengths of 3 x 
18.5 m, and 1 x 10.2 m. Each string consists of 12 pairs of OMs. All OMs face downward, 
with the exception of the OMs of the second and eleventh layers, which look upward. The 
distance between downward oriented layers is 6.25 m, the distance between layers facing 
to each other (layers 1/2 and 10/11) is 7.5 m, the distance between back-to-back layers 
(2/3 and 11/12) is 5.0 m. 

Results presented below are from 70.3 days data taking with NT-96, starting at April 
13, 1996. During the first month, 47 out of all 48 channels have been operating. In the 
following, nine channels failed, so that NT-96 consisted of only 38 operating channels at 
the end of this period. 

3 Identification of Upward Muons by Full Recon- 
struction 

3.1 The Reconstruction Procedure 

The reconstruction algorithm is based on the assumption that the light radiated by the 
muons is emitted exactly under the Cherenkov angle (42 degrees) with respect to the 
muon path. The model of single minimum ionizing muon is a simplification, since the 
direction of shower particles accompanying the muons is smeared around the muon direc- 
tion. Atmospheric muons may occur in bundles. Also, light may be scattered on its way. 
In Lake Baikal, however, the latter effect is negligible since the effective scattering length 
is much larger than the absorption length [|IJ. The reconstruction procedure consists of 
the following steps: 

1. A first quality analysis of the event which a) excludes events far from being described 
by the model of a minimum ionizing muon, and b) finds a first guess for the \ 2 
minimization. 

2. Determination of the muon trajectory based on the minimization of the function 

N hit 

Xt = J2 fa u o, v , *o) - U) 2 /oli 

8=1 
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Here, U are the measured times and Tj the times expected for a given set of track 
parameters. Nh it is the number of hit channels, a t i are the timing errors. All a ti 
are set to 5 nsec in the following procedures. A set of parameters defining a straight 
track is given by 6 and <fi - zenith and azimuth angle of the track, respectively, uq 
and v o - the two coordinates of the track point closest to the center of the detector, 
and t ~~ the time the muon passes this point. For a fit allowing a three-dimensional 
reconstruction, >6 hits on > 3 strings are necessary (trigger 6/3). 

3. Rejection of most of the bad reconstructed events with the help of final quality 
criteria. 

In the initial quality analysis (step 1), an event has to pass the following criteria: 

a) The time difference Ai^ for two channels % and j must obey the following condition: 

| Atij | cos?7 < Tij/c + 5, (1) 

with i] being the Cherenkov angle (42°) and the distance between channels. 5 = 
5 nsec accounts for the timing error. 

b) For any two channels % and j on the same string, a region of zenith angles 6 mm — 9 max 
is determined which is allowed by the observed time differences Atif 

r At - 

cos (6 mm + v )< cosr]^-^ < cos(9 max - rj) (2) 

Zj z% 

Here z^y Zj c^re the z coordinates of the channels. If the regions of possible zenith 
angles for all pairs along a string do not overlap, the event is excluded. 

c) Assuming the event is caused by a naked muon, for every channel one can de- 
fine a range of distances to the muon which is consistent with the measured am- 
plitude of the channel. From this, for every channel pair the minimal (At mm ) 
and maximal (At max ) allowed time difference are defined, in dependence on the 
distance between the channels and the amplitudes. If for any pair the condition 

At min < Af exp < ^max j g violated) the event jg re j eC ted. 

Seventy percent of the events selected for reconstruction (trigger 6/3) pass these criteria 
in both the experimental and the MC sample. In the case of MC generated neutrino 
induced muons, the rate is larger (80%) due to the absence of muon bundles. 

After the minimization (step 2) we apply final quality cuts (step 3). Apart from the 
traditional cut on Xt > cu ^ s 011 the following parameters are applied: 
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• the probability of non- fired channels not to be hit, Pnohit, and fired channels to be 
hit, P hit . 

n m 

Phu = {UPi) 1/n Pnom=(T[Pi) 1/m ( 3 ) 

i=l j=l 

where pi is the hit probability for the ith of all n fired channels, and pj the nohit 
probability for the jth of all m non-fired channels. 

• the correlation function A corr of measured amplitudes to the amplitudes expected: 

n 

A corr = J2( A i~ A)(ai- a)/a a a A (4) 

i=l 

with Ai being the amplitude expected for channel i on the base on the fitted track 
hypothesis and a minimum ionizing muon, tij the measured amplitude and a a and 

a a the corresponding dispersions. 

• an amplitude x\ defined similar to the time Xt defined above. 

Figs.0 and |] show as an example the distributions of Put • Pnohit and the amplitude correla- 
tion Aco rr - for downward muon Monte Carlo events compared to the total reconstructed 
experimental sample, and for fake events compared to true upward muons (both MC). At- 
mospheric muons have been generated according to [[H]] (see for details |lj where also the 
Baikal detector simulation is described), for atmospheric neutrinos the Volkova spectrum 



1 7 has been used. 



To guarantee a minimum lever arm for track fitting, we reject events with a projection 
of the most distant channels to the track (zdist) smaller than a certain length. Due to 
the small transversal dimensions of NT-96, this cut excludes zenith angles close to the 
horizon, i.e., the effective area of the detector with respect to atmospheric neutrinos is 
decreased considerably. 

We have chosen the following cuts to separate well reconstructed tracks: 
1: Xt/NDF < 3.0 

2: Phit ' Pnohit > 0.15 

3: xl/NDF < 2.0 

4: amplitude correlation A corr > 0.1 

5: z dist > 35 m 

Fig.|] shows the effective area for reconstructed upward muons with > 9 hits on >3 strings, 
and for those passing additionally cut 1-4 and 1-5, respectively, as function of the zenith 
angle. 
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Figure 2: Distribution of P^u • Pnohit- Left: downward muon Monte Carlo events (dashed) and 
reconstructed experimental events (full line). Right: MC fake events (dashed) and true MC 
upward muons (full line). 
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Fi gure 3: Distribution of the amplitude correlation A corr . Left: downward muon Monte Carlo 
events (dashed) and reconstructed experimental events (full line). Right: MC fake events 
(dashed) and true MC upward muons (full line). 
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3.2 Results Obtained with the Standard Analysis 

The efficiency of all criteria was tested using MC generated atmospheric muons and up- 
ward muons due to atmospheric neutrinos. 1.8 ■ 10 6 events from atmospheric muon events 
(trigger 6/3) have been simulated, with none of them passing all cuts and being recon- 
structed as upward going muon. Table 1 shows the number of events surviving various 
trigger/cut levels (column 1) - column 2 for the experimental data, column 3 for the 
Monte Carlo sample of downward muons, column 4 for the same sample normalized to 
70.3 days, and column 5 for upward muons from atmospheric neutrinos, again normalized 
to 70.3 days. 

Rejection of all events with less than 9 hits results in a small decrease of the neutrino 
sensitivity (see table 2), but reduces the background by a factor of 4. This corresponds 
to the lowest curve in fig.|4]. 

We have reconstructed 2.0 • 10 7 events taken with NT-96 from April to August 1996 (70.3 
days). Nine events were reconstructed as upward going muons, passed all quality cuts 
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and triggered at least 9 channels at 3 strings. This compares to 8.7 events expected from 
atmospheric neutrinos. Fig. 5 displays one of the neutrino candidates. Top right the times 
of the hit channels are shown as as function of the vertical position of the channel. At each 
string we observe the time dependence, characteristically for upward moving particles. 



Table 1: Number of events surviving various trigger/cut levels 



Trigger /cut 


Experiment 


MC atm n 


MC atm n 
70 days 


MC atm v 
70 days 


6/3 


2.0 ■ 10 7 


1.8 • 10 6 


2.5 • 10 7 


51.6 


(l)+upward rec'd 


7.2 • 10 5 


5.5- 10 4 


7.7- 10 5 


49.7 


(2)+quality cuts 1-4 


5088 


454 


6361 


17.3 


(3)+quality cut 5 


30 







9.9 


(4)+9/3 


9 







8.7 



Table 2: The fraction of events after the quality cuts relative to the number of events 
surviving pre-criteria and reconstruction for the triggers 6/3 and 9/3 



Trigger cond. 


Experiment 


MC atm /i 


MC fi from v 


6/3 


0.19 


0.21 


0.20 


9/3 


0.044 


0.056 


0.175 



Applying eq.0 not only to pairs at the same string, but to all pairs of hit channels, one can 
construct an allowed region in both 9 and 0. For neutrino events this region is requested 
to include regions below the horizon. This is demonstrated at the bottom right picture of 
fig.5. The same holds for the other eight events, one of which is shown in fig.|6|a). Fig.^), 
in contrast, shows an ambiguous event not passing all quality cuts and giving, apart from 
the upward solution, also a downward solution. In this case we assign the event to the 
downward sample. 

The resulting angular distribution is presented in fig.[7|a). Fig.[7|b) shows the measured 
angular distribution of the 9 neutrino candidates weighted with the acceptance and com- 
pares it to the spectrum of atmospheric neutrinos which was taken from ||17|| . Table 3 
gives the number of hit channels and of hit strings, and the zenith angle for the nine 
neutrino candidates. The average number 12 of hit channels agrees within errors with the 
expected number of 11.2, the mean Xt is 0.91, only one event has a \\ value (2.45) close 
to the cut limit of 3. 
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Figure 5: A "gold plated" 19-hit neutrino event. Left: Event display. Hit channels are in 
black. The thick line gives the reconstructed muon path, thin lines pointing to the channels 
mark the path of the Cherenkov photons as given by the fit to the measured times. The sizes 
of the ellipses are proportional to the recorded amplitudes. Top right: Hit times versus vertical 
channel positions. Bottom right: The allowed 9/<p regions (see text). 
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Figure 6: left: An unambiguously reconstructed 10-hit neutrino candidate, right: An ambiguous 
event reconstructed as a neutrino event (dashed line) but not passing all quality cuts. This event 
has a second solution above the horizon (solid line) and was assigned to the sample of downward- 
going muons. 



Table 3: Parameters for the nine upward muon candidates 



# 


N hit 


A' string 


9 


1. 


14 


4 


134.7° 


2. 


19 


4 


152.3° 


3. 


10 


3 


120.4° 


4- 


16 


4 


145.5° 


5. 


11 


4 


121.8° 


6. 


9 


4 


138.6° 


7. 


9 


3 


132.6° 


8. 


11 


4 


136.9° 


9. 


9 


4 


125.3° 



11 




cos 8 cos 



Figure 7: left: Experimental angular distribution of events satisfying trigger 9/3, all final 
quality cuts and the limit on Zdist ( see text), right: The angular distribution below cos# = 
-0.4, corrected for the acceptance. The full line is obtained from the Volkova spectrum for 
atmospheric neutrinos [fL7f . 



4 Identification of Nearly Vertically Upward Moving 
Muons 

Different to the standard analysis, the method presented in this section relies on the 
application of a series of cuts which are tailored to the response of the telescope to nearly 
vertically upward moving muons. The cuts remove events far away from the opposite 
zenith including fake events. Fake events are mostly due misinterpreted atmospheric 
muons with true zenith angles close to the horizon, and to pair and bremsstrahlung 
showers below the array which may give rise to upward traveling light fronts. The method 
does not aim a full reconstruction of all spatial parameters. Therefore, also events with 
hits at only one or two strings are included in the analysis. The candidates identified by 
the cuts are afterwards fitted in order to determine the zenith angle. 

4.1 Definition of the cuts 

In order to optimize the cuts with respect to the signal-to-noise ratio, the following Monte 
Carlo samples have been generated and analyzed: 

• Atmospheric neutrinos below horizon (4.4 • 10 3 triggered events corresponding to 53 
years data taking). 
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• Minimum ionizing atmospheric muons with zenith angles > 60° (1.5 • 10 6 triggered 
events corresponding to 3.6 years data taking). 

• showers generated by downward muons close to the array. Showers have been gen- 
erated in a 160 x 160 x 100m 3 box which includes the detector center. The energies 
ranged from 10 GeV to 10 TeV and were normalized to atmospheric muons (energy 
and angular distribution, and production probability for bremsstrahlung, 3.7 • 10 7 
triggered events corresponding to 10 years data taking) 

We included all events with >4 hits along at least one of the four strings. 
4.1.1 Timing 

The first criterion cuts on the times measured along each of those individual strings with 
more than 1 hit. It applies to all combinations of hit channels along the string: 

cut 1: |(ij — tj) — Azij/c\ < a ■ Azij + 25, (i < j). 

The ti,tj are the arrival times at channels and Az^ is their vertical distance. With 

<5 = 5nsec accounting for the timing error, the condition |(ij — tj) — Zy/c\ < 25 (i.e. a=0) 
would cut for a signal traveling vertically upward with the speed of light, c. 

Setting a to 1 nsec/m, the acceptance cone around the opposite zenith is slightly increased. 
Fig.|3| shows the effective area for muons from atmospheric neutrinos as a function of cos 9 
after cut 1. All following cuts are applied to the sample of events passing cut 1. 



S 




-1 -0.8 -0.6 -0.4 -0.2 

cos 

Figure 8: Effective area for muons from atmospheric neutrinos passing cut 1. 
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4.1.2 Event Length 



Events close to the vertical tend to be much more elongated than hit patterns due to 
showers close to the array, and also than fake downward muons, see fig.||]. The following 
cut is the next powerful after cut 1 with respect to fake event rejection: 



cut 2: 



J vert 



\ibot — hop + 1| > 



with i top and % ot being the numbers of the highest and the lowest floor, respectively, which 
contain hit channels. Numbers are counted from top to bottom. This cut defines a rather 
sharp energy threshold at E„ ps 10 GeV. 



V 

1 10' 



t io 3 

3 



io- 



10 1 



10 



showers 



_r 



signal muons 



background muons 



4 5 6 7 



10 11 12 

L vprt 



Figure 9: Distribution of the event length L ev for events from atmospheric neutrinos (full line) 
as well as fake events from naked downward muons (dashed line) and showers (dotted line), 
respectively. 



4.1.3 Maximum Amplitude 



The distribution of the maximum amplitude in an event, A max , differs strongly between 
atmospheric muons (mean value 23 photoelectrons), naked- muon fake events and shower 
events. 



cut 3: 



< 50 photoelectrons, 

accepts 90 % of the signal and removes about 50 % of the shower events (for the parameter 



distributions underlying this and the following cut values, see [11 
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4.1.4 Center of Gravity 



Naturally, events due to showers which emit light traveling upward within the detector 
deposit most light in the lower part of the array. From the distribution of the center of 
gravity of hit channels for signal and for shower events (not shown), we derive a cut on 
the z coordinate of the center of gravity: 

cut 4: COG z = E,=i(A • z^/ElMi) > 20 m, 

with Ai being the amplitude of the ith channel and N the number of all hit channels. 

4.1.5 Causality Condition 

The timing cut 1 refers to the time differences between hit channels along one string. 
However, one can also relate the hits at different strings to each other assuming the 
particle is moving nearly straight upward. For this reason, we define a parameter 

t — + 4- tt0p — ^ bot (v — y \ 

l basic,i l bot < V i ^bot) 

%top %bot 

Here, Zbot and Zt op are the z coordinates of the lowest and the highest fired channel at a 
string with Nhu > 4, respectively, and Uot and t top are the corresponding times. Zi is the 
z coordinate of a hit channel at another string. In the case of a straight upward muon, 
the time U measured for this channel should be close to tbasic,i- 

From the distribution of the largest values of \U — U MS ic,i\ found for each event, we derive 
the following "causality cut": 

cut 5: t — max(|tj — tbasic,i\) < 60 nsec. 

4.1.6 Minimum time difference between top and bottom channels 

According to cut 2, (L vert > 8), the distance between top and bottom hit channels should 
be at least 50 m. This corresponds to a travel time of 166 nsec for a relativistic muon. 
However, the actual difference in arrival times may be slightly lower due to light scattering, 
shower particles accompanying the muon and experimental errors. We request that for 
all combinations of top and bottom hits from various strings 

cut 6: t tot = wm(t top>i - t botJ ) > 150 nsec 

where denotes the string number. 
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4.2 Application to Experimental Data 



Within 70 days of effective data taking, 8.4 • 10 events with the muon trigger Nha > 4 
have been taken. 8608 of them have >4 hits at one string and pass cut 1. Table 4 
summarizes the number of events from all 3 events samples (MC signal and background, 
and experiment) which survive the subsequent cuts. The experimental sample consists 
of 8.4 • 10 7 events taken within 70.3 days effective live time with the hardware trigger 
Nhu > 4. The Monte Carlo samples are normalized to the same time interval. 

Table 4: The expected number of atmospheric neutrino events and background events, 
and the observed number of events after cuts 1-6. 



after cut N° -> 


1 


2 


3 


4 


5 


6 


atm. is, MC 


11.2 


5.5 


4.9 


4.1 


4.8 


3.5 


background, MC 


7106 


56 


41 


16 


1.1 


0.2 


experiment 


8608 


87 


66 


28 


5 


4 



The disagreement between background Monte Carlo and experiment in the first columns 
is mostly due to the fact that the Monte-Carlo was simplified in order to save CPU time 
and does not include low energy showers (E s h ower < 10 GeV) accompanying the muon. 
Whereas the number of MC background events continues to fall with additional cuts, 
experimental data and MC signal events apparently reach a kind of asymptotic value 
which is only slightly affected by the last cut. The parameters of the four experimentally 
selected events are given in table 5 together with the results of the zenith angle fit. We 
note that all parameters values fall in the most probable range for signal events which 
had passed the time cut. Moreover, they are not close to the cut limits, i.e., a reasonable 
variation of the cuts would not change the number of events passing the cuts. We conclude 
that we have separated 4 atmospheric neutrino candidates with expectation values of 3.5 
for the signal itself and 0.2 for the background from fake events. 

Table 5: Fit results for the four upward muon candidates surviving all cuts 1-6 



# 


N hit 


■N string 


Lev 


A 

(p.e.) 


COG z 
(m) 


t 

(nsec) 


ttot 

(nsec) 


e 

(degree) 


X7NDF 




19 


4 


12 


12.5 


34.8 





190 


152 


0.44 




9 


2 


12 


21.0 


23.5 


38 


270 


161 


0.05 




9 


2 


10 


13.0 


26.0 


22 


202 


162 


0.14 




7 


2 


10 


4.0 


32.0 


23 


188 


170 


0.39 
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The first of the four events is the gold plated 19-hit event which was also found in the 



full-reconstruction analysis. FigfTCj shows the second event with 9 hits at 2 strings. Fig.|Tl 



gives the angular distribution of the four events compared to the effective area after all 
cuts. With a nearly flat distribution for atmospheric neutrinos over the small cos 9 range 
considered, the observed events should be distributed according to the effective area. In 
fact, the distribution of the four observed events is compatible with this assumption. 




Figure 10: The third event of table 4, with 9 hits at 2 strings and 161° reconstructed zenith 
angle. 
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Figure 11: The effective area of NT-96 after cuts 1-6 (histogram) compared to the reconstructed 
angles of the four events passing these cuts. 



4.3 Limit on the Flux of Muons from the Center of the Earth 

One class of candidate particles for cold dark matter has an interaction strength typically 
for the weak scale. Weakly Interacting Massive Particles (WIMPs) would have been 
produced during the big bang and are diffusing through the universe. They are expected 
to be scattered off coherently by nuclei, e.g. in the Sun or the Earth. Loosing energy they 
eventually fall below the escape velocity and are gravitationally trapped at the center 
of these bodies. Once accumulated, annihilation processes set on. The WIMP density 
builds up until equilibrium is reached. The leptonic decays due to annihilation into heavy 
quarks or W-pairs yield neutrinos of energies in the GeV-to-TeV range, depending on 
the WIMP mass. The most favoured candidate for WIMPs is the neutralino. In the 
Minimal Supersymmetric Standard Model (MSSM), neutralinos are linear combinations 
of the neutral gauginos B, Ws, and of the neutral higgsinos H°, H^. The lightest of these 
combinations would be the stable candidate for cold dark matter. The neutralino capture 
rate in celestial bodies, the annihilation cross sections, the fragmentation functions into 
fermions, gauge bosons and Higgs particles depend on the various parameters of the 
MSSM (see [112, [T3j and references therein). 



Underground detectors like Baksan, MACRO and Kamiokande set upper limits on the 
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flux of muons from the center of the Earth of a few times 10 14 cm 2 sec 1 M, EL6L [15 



These limits already exclude a certain region of the SMMS parameter space. Neutralino 
masses below 25 GeV are excluded by LEP. 

Table 5 shows the number of events detected in NT-96 within four zenith angle regions 
about the opposite zenith, and compares it to expectations from atmospheric neutrinos. 
From the non-observation of a significant excess, the flux limits given in the third column 
are obtained. They are compared to the limits given by Baksan, MACRO and Kamiokande 
in the last 3 columns. In spite of the effective area of NT-96 for nearly vertically muons 
being about 1000 m 2 (see fig. [TT|) , the limits obtained are 4-7 times wore than the best 
published values, since the only 70 days have been analyzed until now. However, the 
result illustrates the potential of underwater experiments with respect to the search for 
muons due to neutralino annihilation. 



Table 5: The number of events detected ("Data") and expected from atmospheric neu- 
trinos ("Bg") for NT-96, as well as the 90% C.L. upper limits on the muon flux from the 
center of the Earth for four regions of zenith angles obtained NT-96, Baksan, MACRO 
and Kamiokande. 



Zenith 
angles 


Data 


Bg 


Flux limit (10 14 • cm 2 sec 1 ) 


NT-96 
> WGeV 


Baksan 
> IGeV 


MACRO 
> 1.5GeV 


Kamiokande 
> 3GeV 


> 150° 


4 


3.7 


11.0 


2.1 


2.67 


4.0 


> 155° 


3 


2.6 


9.3 


3.2 


2.14 


4.8 


> 160° 


2-3 


2.3 


5.9-7.7 


2.4 


1.72 


3.4 


> 165° 


1 


1.3 


4.8 


1.6 


1.44 


3.3 



5 Conclusions 

Twelve neutrino candidates have been separated from 70 days effective lifetime of the 
four-string detector NT-96 in Lake Baikal. Nine of them have hits at > 3 strings and 
are fully reconstructed (zenith and azimuth angle). The other three events have hits at 
2 strings but still yield an unambiguous zenith reconstruction. The observed number as 
well as the angular distribution coincide with expectations from atmospheric neutrinos. 

This is the first successful experimentum crucis for the operation of underwater neutrino 
telescopes. From the agreement with expectations, we derive an upper flux limit of 
1.1 ■ 10~ 13 cm -2 sec -1 for muons with zenith angle > 150° and >10GeV in excess 
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of those expected from atmospheric neutrinos. Such muons might have been produced in 
neutralino annihilations in the center of the Earth. 

At present, the NT- 200 telescope is taking data. It contains twice the number of optical 
modules as NT-96. From the fully operational NT-200 array, about one neutrino event 
per day is going to be separated. NT-200 will be primarily used to study atmospheric 
neutrinos and to search for upward muons from WIMP annihilation. 
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